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We show that the whole range of RHIC data for hadron production in d-Au collisions 
is compatible with geometric scaling. To establish the scaling violations expected from 
small- x evolution a larger kinematic range in transverse momentum and rapidity would 
be needed. We point out that the fall-off of the pt distribution of produced hadrons at 
large pt is a sensitive probe of small- a; evolution especially at the LHC. 

It is clearly observed that the small- a; DIS data show the property of geometric scaling 
(GS) [2J, i.e. the cross section depends on the combination Q 2 /Q 2 (x) only, where Q s (x) is 
referred to as the saturation scale. On the other hand, geometric scaling is a feature of the 
asymptotic solutions of nonlinear evolution equations, such as the BK equation [5]. Hence, 
geometric scaling is seen as a strong indication for small- a; gluon saturation. 

A phenomenological study of DIS data using a model for the dipole cross section was 
performed by Golec-Biernat and Wusthoff (GBW) [I]. They found that the inclusive HERA 
data at low x < 0.01 could be described well by a dipole cross section a = o'oA f GBw(''( i x), 
where <j$ ~ 23 mb and the scattering amplitude -/Vqbw is given by 

N GB w(n,x) = 1 - exp (~\r 2 t Q 2 s {x)) . (1) 

The x-dependence of the saturation scale is given by 

K/2 , Qo = 1 GeV (2) 
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where the parameters x n ~ 3 x 10~ 4 and A ~ 0.3 are fitted to the small- a; data. The 
amplitude |T]) depends on x and r t (the transverse size of the dipole) only through the 
combination r 2 Q 2 (x). In DIS, this directly results in a scaling cross section. Hence, in the 
dipole picture, GS is equivalent to a dependence of the amplitude on r\ Q 2 (x) only. 

Despite that GS is expected already sub-asymptotically in a growing range around Q s 
[5], the values of x probed at recent colliders may be not small enough to expect GS and 
it may be more important to test violations of GS to establish small- a; evolution. In order 
to investigate GS violations in the RHIC data, in [U [7] a phenomenological model has been 
put forward. We will refer to this model as the DHJ model. It offers a good description of 
the pt distribution of hadrons produced in d-Au collisions at RHIC in the forward region 0. 

According to Refs. [6l [7] the cross section of hadron production in high-energy nucleon- 
nuclcus collisions can be described in terms of the dipole scattering amplitude, 
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a As it turned out the study of Ref. [7] contained an artificial upper limit on the x\ integration to exclude 
large X2- Without this cut, the larger pt data for = 0, 1 are in fact not well-described by the DHJ model. 
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Here Np describes a quark scattering off the nucleus, while Na applies to a gluon. The 
parton distribution functions f q / p and the fragmentation functions D^/q are considered at 
the scale Q 2 = p 2 , which we will always take to be larger than 1 GeV 2 . The momentum 
fraction of the target partons equals x 2 = £1 exp(— 2?//j). We can to good approximation 
neglect finite mass effects, i.e. we equate the pseudorapidity 77 and the rapidity yh and use 
xf = V p 2 + m 2 /^fs exp(?7) f=a pt / \fs exp(yh) ■ Finally, there is an overall yh dependent 
if -factor that effectively accounts for NLO corrections. It should be noted that the scaling 
properties of the dipole scattering amplitude are not directly visible in the hadron production 
data, due to its convolution with the parton distribution and fragmentation functions. 
The dipole scattering amplitude of the DHJ model is given by [5J [7] : 

N A (q u x 2 ) = J d 2 r t [l - exp (-^Q 2 ^)) 1 ^^)] ■ (4) 

Note that 7 is a function of q t rather than r t . This allows one to compute the Fourier 
transform more easily. The corresponding expression Np for quarks is obtained from N A 
by the replacement {r 2 Q 2 s )~< -» {{C F /C A )r 2 Q 2 s )' 1 , with C F /C A = 4/9. The exponent 7 is 
usually referred to as the "anomalous dimension", although the connection between N a /f 
and the gluon distribution inside the nucleus cannot always be made. 
The anomalous dimension of the DHJ model is parameterized as 

7D HJ (gt,S2 )=7. + (l-7.) - —L-J±±L > K 5 

Xy + d^/y+ \ log{crf/Qj(x2))\ 

where y = logl/a;2 is minus the rapidity of the target parton. The saturation scale Q s (x2) 
and A are taken from the GBW model ©. Here Q s includes a larger value of Qo ~ 1-63 GeV 
to account for the size of the nucleus. The parameter d = 1.2 was fitted to the data. This 
choice of 7 leads to a geometric scaling solution at q t = Q s where 7 = 7 S = 0.628 and 
incorporates to a certain extent the violation expected from BFKL evolution for larger q t , 
see [5]. However, an analysis of the BK equation suggests that a smaller value of 7 ss 0.44 
[8] may be more appropriate at Q s where the BFKL equation does not apply. 

Not only a constant 7 would lead to GS but also a 7 that depends on q 2 /Q 2 or r 2 Q 2 
only. To check explicitly whether the RHIC data are compatible with GS, we propose a new 
scaling parameterization of 7 that is similar in form to that of the DHJ model, but does not 
have the GS violating behavior nor the logarithmic rise expected from the BFKL (and more 
generally, BK) equation. The parameterization that we adopted reads 

7gs(^ = qt/Qs) =11 + {I- 7i) - ( f" 1 !. • (6) 

(w a — 1) + 

The two free parameters a and b will be fitted to the RHIC data. There are two major 
differences between the chosen parameterization 7qs © and 7dhj ©■ Firstly, 7qs does 
not depend on the rapidity y explicitly. Therefore the resulting dipole scattering amplitude 
respects geometric scaling. Secondly, 7cg approaches the large q t limit of 1 much faster. This 
will lead to different large momentum slopes of the amplitude (© and therefore to different 
predictions for the large p t slope using Eq. © . For large w the exponential function can be 
expanded and the fall-off of the dipole scattering amplitude © is given by [S] , 

— 7 % ._ 2 . for 7 of Eq. © 
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for 7 of Eq. © 
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Empirically, we find that the pt distribution falls off even faster. We note that the fall-off 
with p t is not determined by the size of the scaling violations. In order to observe such 
violations, one has to study both the yh and p t dependence over a significantly large range. 
Let us mention explicitly that our parameterization is not meant to replace other, physically 
better motivated models but to investigate in a general way which conclusion can really be 
drawn unambiguously from the RHIC data in the central and forward regions. 

In Fig.[T]a) we show our estimate for dNh/ '(dyuCppt) that follows from the integral in Eq. 
([3]) by using r )Gs{w) ([6]). All p t distributions of produced hadrons measured at RHIC in d-Au 
collisions are well described. At the saturation scale we have chosen here for 7gs(' u; = 1) = 7i 
the same value 7 S = 0.628 as in the DHJ model. We also take the same parameterization 
of Q s (x). We obtain the best fit of the data for a = 2.82 and b = 168. As mentioned, this 
LO analysis requires the inclusion of a K-i&ctov to account for NLO corrections, which are 
expected to become more relevant towards central rapidity. The p t -independent K factors 
we obtain for yh = 0, 1, 2.2, 3.2, 4 are for our model equal to K = 3.4, 2.9, 2.0, 1.6, 0.7 and 
for the DHJ model K = 4.3, 3.3, 2.3, 1.7, 0.7. For further details of the calculation see [9j. 




5678v 02468 
Pt[GeV\ w = „/(J, 



Figure 1: a) Transverse momentum distribution of produced hadrons in d-Au collisions as mea- 
sured at RHIC (symbols) for various rapidities yh- To make the plot clearer, the data and the 
curves for yh = 0,1 and 2.2 are multiplied with arbitrary factors, namely 16, 4 and 2, respectively, 
b) Various fits of 7gs(w), which describe the RHIC data equally well. For comparison we show 
curves representing yunj(w, y(w, yh)) at different rapidities yh- 

From this analysis we can conclude that a GS dipole amplitude is completely compatible 
with the data and therefore the conclusion that GS violations are observed at RHIC cannot 
be drawn. Of course, a scaling violating amplitude, i.e. a 7 that depends on w and the 
rapidity y explicitly, is not ruled out by the data either. What can be concluded further is 
that the logarithmic rise of 7 resulting from the BFKL evolution incorporated in the DHJ 
model is ruled out in the central region, see Fig. Q]a). However, where the DHJ model starts 
to deviate from the data xi becomes larger than 0.01, although Q s then is still larger than 
in DIS at x = 0.01. If one were to exclude the central rapidity RHIC data in the model fit, 
one could also obtain a scaling model with a logarithmically rising, or even constant, 7. 

To indicate how much 7 is constrained by the RHIC data, Fig. [1] b) shows various 
7gs(w)'s that describe the available data equally well. They are all parameterized as in 
Eq. ^ with different a and b values. Furthermore, we added lines representing 7dhj © 
for different rapidities. For this one needs to express the rapidity of the target parton y in 
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terms of w and yh, see [5] for details. It should be mentioned that the region below Q s , 
where the parameterization of 7dhj is not smooth, is hardly probed at RHIC. Clearly, 7 is 
less well determined close to the saturation scale than in the dilute region. This is because 
the integrand entering the dipole scattering amplitude (j4|) is only weakly dependent on 7 
around the saturation scale r = \/Q s . In addition, the forward data (yh = 3.2 and 4) are 
essentially sensitive only to 71, since they probe the region where w is close to 1. Therefore, 
the rise of 7 with w is effectively constrained only by the data for yh = 0, 1. It is apparent 
from Fig.[T]b) that the DHJ model fails for larger w, i.e. larger p t , in the central region but 
not in the forward region where the probed values of w are below 2. 

Where the DHJ model curves deviate from the RHIC data, the probed a^-values are not 
very small. However, at LHC, due to the higher energies, the region of small x 2 extends to 
a much larger range of pt, so that the predictions of the DHJ model and the new one will be 
different even at small x 2 . It has to be mentioned that the calculation of the pt distribution 
using our scaling model should not be seen as a physically motivated prediction. However, a 
comparison of the estimate using 7gs fitted to RHIC data with the estimate from the DHJ 
model allows drawing conclusions anyway. As we will explain, the estimates are so different 
that the LHC should be able to rule out one of these models at much smaller x. Hence, 
the LHC can answer the question why the DHJ model fails in the central region at RHIC. 
Either because the probed x values are not small enough or because some expectations from 
small- x evolution has to be modified. For smaller p t the predictions of the two models are 
expected to be comparable since only the region of small values of w is probed where 7dhj 
and 7gs are similar, see Fig. [T]b). In the very forward region, i.e. yh ps 6 — 8, only this 
region is tested. However, there is quite a large range where the probed values of X2 are 
small but the predictions are clearly different. One expects (see Fig. [T]b)) differences in the 
central region between the two models for values of w > 3, i.e. at RHIC for p t > 2.5 GeV. 
Since Q s is larger at LHC such differences show up at larger momenta around 5 GeV. 
Even in the central region such momenta imply values of X2 smaller than 0.001. Hence, a 
measurement of the slopes at moderate rapidities yh at LHC would allow a discrimination 
between the two models in a region where small- x physics is expected to be applicable. 
The slower fall-off of the pt distribution in the DHJ model is a direct consequence of the 
logarithmic rise of 7 towards 1. Since such a behavior of 7 is a generic signature of BFKL 
evolution, these measurements offer the possibility of testing whether such small- x evolution 
is actually relevant at present-day hadron colliders. Estimates of the pt distributions for 
hadron production at LHC using both models can be found in [5]. 
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